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Overview 


Low  frequency  electrochemical  impedance, 
EIS,  as  a  metric  for  corrosion  protection  by 
coatings 

Modeling  and  the  decisions  made  about 
how  to  do  it. 

Apply  same  modeling  principles  as  we  have 
used  to  model  appearance,  toughness 

Examples 
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Electrochemical  Impedance  Spectroscopy 
as  a  measure  of  coating  performance 


Measure  electrical  impedance  of  a  coating  across 
its  thickness 

-  Between  conductive  electrolyte  on  top  and  metallic 
substrate 

-  Intact  coatings  have  very  high  impedance  (both  resistance 
and  capacitance) 

-  Low  frequency  impedance  (resistive  component) 
deteriorates  upon  weathering  until  it  diminishes  to  the 
value  of  a  bare  metallic  substrate 

-  Both  resistance  and  capacitance  can  be  calculated  with 
effective  media  theory  based  on  the  concentration  of  water 
or  ions  incorporated  into  the  film 
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How  do  water  and  ions  get  through 

a  coating? 


Flux/ area/ (concentration  difference)  is  the 
permeability  coefficient  across  a  film  (membrane) 
(measured  at  a  steady  state) 

-  Permeability  coefficient  =  diffusion  coeff.  x  solubility  coeff. 
So  we  need  to  understand  diffusion  and  solubility 
Solubility 

-  Tend  to  find  only  2  -  4%  water  at  saturation  in  an  intact 
coating  polymer 

-  It  is  very  difficult  to  find  any  chlorine  or  other  ions  in  the 
cross-section  of  an  intact  coating,  even  after  prolonged 
exposure 

-  Therefore  not  much  effect  will  be  seen  from  dissolved 
matter  -  exclude  from  modeling  (today) 


How  much  corrosive  material 
diffuses  through  a  coating? 


Ionic  radii 

-  Sodium  loses  an  electron,  gets  smaller 

-  Chloride  and  sulfates  gain  electrons  and  get  bigger 


• 

'  ■  '  

Water  molecules  are  small  but  do  we  ever  find  one  by  itself? 
Diffusion 


-  Diffusion  coefficient  of  water  into  polymers  tends  to  be  minute, 
m^/s,  takes  many  hours  to  saturate  (even  to  2-4%) 

-  Diffusion  of  NaCl  has  been  measured  and  is  ~  5-6  orders  of  magnitude 
slower  than  water  [Kloppers,  Belluci,  Latanision,  Brennan,  J.  Appl.  Polym.  Sci. 
48,  2197  -  2205  (1993)  -  more  recent  than  Leidheiser  or  Mayne] 

Can  conventional  ions  get  through  an  intact  crosslinked  polymer? 

-  Not  easily:  Chloride  does  not  move  far  in  molecular  dynamics 
simulations  of  crosslinked  or  thermoplastic  polymers  [Eslami,  Miiller- 
Plathe,  J.  Chem.  Phys.  131,  234904  (2009)] 

-  So  it  is  hard  to  imagine  seeing  the  effects  in  an  EIS  test  after  weathering 
if  the  coating  is  intact 
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•  Therefore  we  must  model  coatings 
that  are  not  intact. 

-  Degradation  produces  a  pitted 
and  degraded  structure, 
through  which  water  etc.  find 
it  much  easier  to  pass. 

-  Changes  to  becoming  a  “semi- 
permeable”  membrane  at  some 
stage. 

•  Intact  vs.  degraded  coatings 
issue  also  contributes  to  why 
accelerated  corrosion  tests  are 
not  predictive. 


Surface  SEM  images  of  films  made  from  1.4  PEGDGE/ehitosan  prepolymerization  solutions,  (b)  2.0 
PEGDGE/chitosan  prepolymerization  solutions.  Published  in:  Bryan  D.  MeCloskey;  Hao  Ju;  Benny  D.  Freeman; 
Ind.  Eng.  Chem.  Res.  DOl:  10. 102  l/ie901 197u  Copyright  ©  2009  American  Chemical  Society 
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EIS  Testing  of  weathered  Materials 


Periodic  evaluation  of  weathered  samples 

When  tested,  the  samples  usually  start  dry,  then  we  put 
aqueous  electrolyte  on  them 

-  Measure  the  impedance 

•  Wait  until  impedance  stabilizes,  but  water  gets  in  quickly  (~  30 
minutes)  and  we  record  a  reduced  impedance  as  the  material  is 
degraded. 

Diffusion  of  water  through  intact  coatings  takes  hours 

-  Sometimes  leads  to  unstable  EIS  results 

•  Water  may  get  concentrated  near  substrate  depending  on  cycle 

Model  here  assumes  that  water  etc.  only  gets  into  the 
degraded,  outer,  portion  of  the  coating 

-  In  fact  we  assume  that  the  saturated  (from  test  electrolyte)  outer 
layer  is  a  dead  short  and  we  calculate  the  impedance  of  the 
remaining  thickness. 
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Model  the  depth  of  degradation 

(protection  is  provided  only  by  intact  material) 


Assume  wet, 
porous  layer  is  a 
short  circuit 

Neglected  effects 
(by  comparison) 

-  Water  will  diffuse 
into  the  coating 
from  its  whole 
perimeter 

-  There  will  be  added 
degradation  around 
perimeter 
according  to  Beer’s 
law  of 

photodegradation 

penetration 


Test  Electrolyte 


Zero  resistance, 
no  capacitance 


Intact 

Coating 


Decreased  resistance 
Higher  capacitance 


Substrate 
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Original  Coating 


Impedance  of  One  layer  of  a  Coating 

System 

We  focus  on  calculating  the  impedance  of  the  intact,  interior 
layer 

Interest  ceases  stops  when  corrosion  starts,  i.e.  no  inclusion 
of  double  layer  at  corroding  metal  surface  etc. 

No  modeling  of  inhibition 

Assume  that  a  useful  model  is  a  resistance  and  a  capacitance 
in  parallel 

Low  frequency  impedance  is  dominated  by  resistive 
component 


Assume  thickness  change  due  to  degradation  is  only  variable 


in  resistance: 


Area 
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Calculating  the  reduction  in  intact  thickness: 
Damage  Accumulation  via  Statistics  of  random 
processes  (physics  and  maths) 


Has  been  used  to  predict  gloss,  color,  toughness,  etc. 

Can  be  approached  via  Monte  Carlo  computer  simulations,  but 
we  pursue  an  algebraic  form  here. 

Random  repeated  events  (UV  photons,  sand  etc.)  produce  a 
damaged  surface  that  becomes  random  after  many  events 

-  thus  we  can  use  a  result  from  Central  Limit  Theorem  (in 
Normal/ Gaussian  distributions)  for  the  standard 
deviation 

Standard  deviation  =  (number  of  events) 

Similarly,  bulk  morphology  will  become  randomly  damaged  if  it 
degrades  by  an  accumulation  of  repeated  events. 
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Application  of  Gaussian  Statistics  to 
Changes  in  Surface  Roughness 


N  (t)  =  Average  number  of  fragments  removed  at  a  point  on  the  coating 
=  Flux  of  photons  x  Absorption  probability 

X  Quantum  efficiency  of  removing  a  length  of  polymer 
X  Time  (t) 

=  (^A.Q.t 
=  constant  x  t 

i.e.  average  thickness  (or  mass)  decrease,  Ad,  is  linear  with  exposure 
time  (often  good  approximation). 

If  ^  =  fragment  size 
Ad  =  N(t)  X  s 
Ad  =  X  £ 

In  this  case,  the  rate  of  degradation  is  not  limited  by  concentration  of  any 

material  -  ‘zero’  order  chemical  kinetics,  but  we  can  use  other  kinetics  (later) 
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f  Thickness  diminishes  and  Surface  Roughness 

increases  with  exposure. 

Roughness 

<j  =  increase  given  by  the  standard  deviation  in  the  average  number  of 
polymer  segments  removed  at  a  location,  after  sufficient  exposure  has 
occurred  to  randomize  the  surface,  (s  =  characteristic  physical  linear 
size  of  fragment) 

Roughness  =  =  s^N  =R^(pAQ.t 


200 


Mean  diminishes  linearly 
with  exposure  time. 


Roughness  increases  as  the 
square  root  of  exposure  time 
in  the  case  of  zero  order 
kinetics. 


1/2 

(Exposure  time) 


/Sd—s.N{t^ 
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Depth  of  degradation  using  “zero  order”  kinetics 

From  an  initially  smooth  surface 


•  Gaussian  (normal) 
distributions 

•  First  pore  through  will  be 
the  determining  event 

-  Not  the  average 

•  Average  thickness 
diminishes  linearly  with 
exposure  time, 

•  But  the  distribution 
spreads  out  as  the  variance 
increases 

•  Calculate  the  intact 
thickness  by  subtracting 
the  deepest  degraded 
thickness  from  original 
value 


Distribution  of  degradation  depth 


Position 

•Effective  thickness  remaining,  depends  on 

-Average  change  (linear  with  exposure) 
-Std.  deviation  (square  root  of  time  ) 
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Substrata 


Options  at  this  stage: 

Degradation  may  be  more  complicated  than 

“zero  order” 

From  before,  reduetion  in  mean  thiekness”  Ad  ,  (f) 

Degradation  Kineties  may  be  autoeatalytie 

-  The  degradation  feeds  on  itself  AQN(t) 

dt 

The  increase  in  damage,  over  the  starting  value 

^Changed  (f)=A?o[exp(^AQi)-l] 


Or,  order 

-  the  supply  of  labile  material  is  finite 


dN{t) 


dt 


(l)A.Q.N{t) 


-  The  increase  in  damage 


^Changed  (0  =  A^0  [l-eXp(-4^ei)] 
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The  same  options  apply  to  the  distribution  in 

degraded  thickness 


If  there  is  a  pre-existing  surface  variance  (thickness  variation),  that  is 

also  Gaussian 

=  (t) 

Zero  Order  (see  prior  example): 

Autocatalytic  (Malthusian)  kinetics: 

{t)=<yLa,  -  l) 

First  order  kinetics 

So  the  behavior  of  the  average  depth  of  degradation  and  the  greatest 
depth  of  penetration  can  follow  other  kinetics 


B.  R.  Hinderliter,  S.  G.  Croll,  “Predicting  Coating  Failure  using  the  Central  Limit  Theorem 
and  Physical  Modelling”  Accepted  for  publication  in  ECS  Transactions 
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Example  of  Modeling 


If  water  gets  in  via  damage  due 
to  ‘zero  order’  degradation 
kinetics 

-  there  will  be  term  in  )  and 
linear  t 

-  As  the  low  frequency  impedance 
(resistance)  is  reduced  (data  from 
Gordon  Bierwagen) 

Thickness  of  the  film,  d,  is 
reduced  by: 

d{t)  —  d^—  {A.t  +  B.4t ) 

Coating  life  is  limited  at 
downturn 

Failure  is  when  water  arrives  at 
substrate 

-  vs.  size  of  thin  spot,  plugging 


Logio  (Impedance, 
ohm.cr^,  0.013  Hz) 


Hinderliter,  B.  R.,  S.  G.  Croll,  et  al. 
(2006)  Electrochimica  Acta  51(21): 
4505-4515. 


N.B.  looks  like  a 
diffusion  process 


Other  aspects  &  details  are  possible  via  this 

approach 


Another  example  of  fit  to  linear  plus 
square  root  time. 


Hinderliter,  B.  R.,  S.  G.  Croll,  et  al.  (2006) 
Electrochimica  Acta  51(21):  4505—4515. 


Monte  Carlo  computer  simulations 

-  Percolating  water  channels  seem  to 
follow  a  log-normal  arrival  time 
distribution 

Effect  of  water  content  in  intact 
layer  (changes  resistivity  and 
dielectric  properties) 

-  Via  effective  media  theories 
Not  big  effects 

Effect  of  form  of  water  inclusions 

-  O.  A.  Stafford,  B.  R.  Hinderliter  and  S.  G. 
Croll,  Electrochimica  Acta,  51(17),  3558- 
3565  (2006) 

Examine  whole  frequency 
spectrum,  i.e.  include  capacitance 
(actually  does  not  add  much) 

Other  chemical  kinetics 

-  Data  is  not  discriminating  enough  to 
choose 
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Summary 


•  Examine  low  frequency  EIS  impedance  as  an  indicator  of  remaining 
corrosion  protection 

•  Assumption  is  that  reduction  in  impedance  is  due  to  increase  in 
surface  porosity  effectively  diminishing  the  intact  thickness 

-  Aggressive  species  do  not  easily  get  through  intact  coating  during  typical 
testing 

•  Properties  of  the  intact  material  do  not  change 

•  Simple  (zero  order)  degradation  kinetics  give  dependence  of  coating 
resistance  that  diminishes  as  exposure  time  and  (exposure  time)h2 

-  Seems  to  describe  much  of  the  behavior 

Future 

-  Models  of  corrosion  protection  (EIS)  can  be  linked  to  other 
coating  properties 

-  Provide  starting  parameters:  A,  Q,  ^and  kinetics 
(chemistry) 
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Thanks 
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Questions? 
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Diffusion  Coefficient  Variations  with  (intact) 

Coating  Thickness 


•  Diffusion  coefficient  and 
span  of  diffusion 
coefficients  increase 
with  thickness 


20  30  40  50  60  70 

film  thickness  (^m) 


- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

20  30  40  50  60  70  80 

thickness  (^m) 


Water  Volume  Fraction  (WVF) 
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20  30  40  50  60  70  80 

thickness  (urn) 


90 


10  20  30  40  50 

film  thickness  (um) 


60 


70 


Topcoat’s  WVF 
increase  slightly 

Primer’s  WVF 
increases 


Simulation  of  EIS  with  pore  structure 


Max:  6.00  Max 


5.8 


‘  5.6 


•  •  5.4 


■  5.2 


4.8 


4.6 


-  4.4 


4.2 


:  0.800 
0.8 


0.6 


0.4 


0.2 


Min:  4.00  Min:  0.200 


time  (min) 

Simulated  SF  EIS  fit  to  nonlinear 
equation  just  like  measurements 


•  Fitting  of  simulated  materials 
response  gives: 

-  Relative  dielectric 

•  Saturated  and  dry 

-  Diffusion  coefficient  (effective) 

-  Water  volume  fraction  based  on 
Brasher-Kingsbury  model 


Contour  lines  indicate  fraction  of  concentration 
and  shading  indicates  relative  dielectric  (based  on 
Hartshorn  1905  EMT) 
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Water  Volume  Fraction 


Finite  Element  Simulations  2D  Pore 
fit  to  B-K  single  frequency  impedance 


•Both  dry  and  saturated  relative  dielectric 
increase  at  high  area  fractions 

•Larger  increase  for  deeper  pores 
•Diffusion  coefficient  increase  as  pore 
density  increases 

Consistent  with  measurements  but  not 
definitive 

1  mm  pore  in  EIS  cell  has  an  are^g 
fraction  of  10’^ 


Pore  Area  Fraction 


EIS  Response  to  composites  shapes  (Finite  Element  Geometry  for 

cylindrical  water  inclusions  (3%)) 

Stafford,  O.  A.,  B.  R.  Hinderliter,  et  al.  (2006).  "Electrochemical  impedance  spectroscopy  response  of  water  uptake  in 
organic  coatings  by  finite  element  methods. "  Electrochimica  Acta  52:  1339-1348. 


Height/radius=l/2 


Height/radius=5 


1e-5 
X  1e-4 


Height/radius=10 
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C/Cno  water 


Capacitance  as  a  Function  of  Water  Volume 
Fraction  for  Spheres  and  Cylinders 

Stafford,  O.  A.,  B.  R.  Hinderliter,  et  al.  (2006).  "Electrochemical  impedance  spectroscopy 
response  of  water  uptake  in  organic  coatings  by  finite  element  methods. "  Electrochimica  Acta 

52:  1339-1348. 


Volume  Fraction  of  Water,  % 
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Slope  determines  lifetime 

•  Slope  depends  on 

-  Weathering  conditions 

•  Accelerated 

•  Natural  weathering 

-  Coating  conditions  (extrinsic  properties) 

•  Surface  treatment/ preparation 

•  Coating  thickness 

-  Coating  properties  (intrinsic  properties) 

•  Binder  crosslink  density,  water  volume  fraction  at 
saturation,  diffusion  coefficient,  PVC 
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Gaussian  (Normal)  Distribution 


•  Mean  =  mode  =  median  =  ju 

•  Variance  =  square  of  standard 
deviation  =  a 


•  Central  Limit  theorem  says  a  ~ 
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Simple  calculation 

•  Epon  828  (n  =  1.2)  (molecular  weight 
350  -  400)  and  Epicure  3115 

-  Molecular  weight  between  crosslinks  = 
930  g/mol. 


H2C 

0 
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Monte  Carlo  Percolation 

The  time  to  percolate  water  to  the  substrate  from  the  surface  is 

a  log-normal  distribution. 


Distribution  of  iog(time  to  substrate) 

8 

- ' - r 

99.5- 

98- 


8 

Bin 
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IZIq  012Hz  Time  under  QUV/Prohesion 
Exposure  (expanded  time) 


a 

N 


1E+10 

1E+09 

1E+08 

1E+07 

1E+06 

1E+05 

1E+04 

1E+03 

1E+02 

0  5  10  15  20  25  30  35  40  45  50  55  60  65  70  75  80  85  90  95  100105110115120125130135140145150 

30 


n 


-Yellow  Primer 

Self-Priming  Topcoat 

-  Primer  +  Flat  Topcoat 

Primer  +  Glossy  Topcoat 

Fit  Primer  +  Glossy 
Topcoat 


Time  (Week) 


